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ABSTRACT 

We describe Chandra High-Energy Transmission Grating Spectrometer observations of RT Cru, 
the first of a new sub-class of symbiotic stars that appear to contain white dwarfs (WDs) capable 
of producing hard X-ray emission out to greater than 50 keV. The production of such hard X-ray 
emission from the objects in this sub-class (which also includes CD —57 3057, T CrB, and CH Cyg) 
challenges our understanding of accreting WDs. We find that the 0.3 - 8.0 keV X-ray spectrum of RT 
Cru emanates from an isobaric cooling flow, as in the optically thin accretion-disk boundary layers 
of some dwarf novae. The parameters of the spectral fit confirm that the compact accretor is a WD, 
and they are consistent with the WD being massive. We detect rapid, stochastic variability from the 
X-ray emission below 4 keV. The combination of flickering variability and a cooling-flow spectrum 
indicates that RT Cru is likely powered by accretion through a disk. Whereas the cataclysmic variable 
stars with the hardest X-ray emission are typically magnetic accretors with X-ray flux modulated at 
the WD spin period, we find that the X-ray emission from RT Cru is not pulsed. RT Cru therefore 
shows no evidence for magnetically channeled accretion, consistent with our interpretation that the 
Chandra spectrum arises from an accretion-disk boundary layer. 
Subject headings: binary stars: general — white dwarf: accretion - X-rays 



1. INTRODUCTION 

Symbiotic stars are interacting binaries in which a hot, 
compact star accretes from the wind of a red-giant com- 
panion. Although a few symbiotics contain neutron- 
star accretors (e.g. GX 1+4, 4U 1700+24 4U 1954+31 
IGR J16194-2810: iDavidsen et all Il977t iMasetti etall 
120021 iGallowav et all 120021 : IMasetti et al.ll2007aU bD. the 
accreting compact object is usually a white dwarf (WD). 
Typical binary separations are on the order of AU, with 
orbital perio ds on the order of a few hundred da ys to a 
few decades (|Kenvodll986tlBelczvhski et al.ll2000T ). Sym- 
biotics can thus be thought of as very large cousins of 
cataclysmic variables (CVs). The accretion rate onto the 
WD appears to be high enough in most symbiotic sys- 
tems that accreted material is burned quasi-steadily in a 
shel l on the WD surface, p roducing a high UV luminos- 
ity ((Sotoioikriran [111 [OrioiLan[2003). Although 
accretion disks a re likely to exist around the WDs i n 
some symbiotics (|Liviolll997HSokoloski fc Kenvonll2003T ). 
there is little direct evidence for these disks. Finally, the 
red-giant wind produces a dense nebula that surrounds 
the system. 

Most symbiotic stars with detectable X-ray emis- 
sion display soft, or supersoft, thermal X-ray spectra. 
As in the supersoft X-ray sources, the lowest-energy 
X-rays could emanate directly fro m material burnin; 
quasi-steadily o n the WD surface ([Jordan et alJ 1199 
lOrio et all 120071). Based on a survey of symbiotics with 
RO S AT, iMuerset et"aH (|1997f) proposed that symbiotics 
be classified according to the hardness of their X-ray 



spectra. They labeled sources with supersoft spectra 
a-types, sources with the slightly harder spectra likely 
to arise from collision of the red-giant and white-dwarf 
winds /3-types, and systems with the hardest spectra 
that might be indicative of neutron stars 7-types. Al- 
though more recent observations using the broader X-ray 
coverage and greater sensitivity of Chandra and XMM- 
Newton have shown that some symbiotics do not fit into 
the si mple a/ 0h class ificat ion scheme (e.g., Z An d and 
o Cetk lSokoloski et al.ll200llKarovska et al.ll2005f ). most 
still appear to produce primarily soft X-rays {E < 3 
keV). 

With the advent of the sensitive hard X-ray detec- 
tors on the Swift and INTEGRAL satellites, a new pic- 
ture has emerged. Some symbiotic stars can produce 
X-ray emission out to greater than 50 keV. Such hard X- 
ray emission has so far been detecte d from 4 symbiotics 
thought to harbor WD s - RT Cru (IChernvakova et al.l 
2005Y IBird et all [200l . T CrB (iTueller et alJ 1200a 



Electronic address: 
Electronic address: 



K.jmluna@astro.iag.usp.br 



jeno@astro.columbia.edu 



Luna et alJl2007f). CH Cve fMukai et al II2007T) and CD 
57 3057 (|Masetti et al.112009 : IBird et al.ll2007| ). Although 
the origin of this hard X-ray emission is not known, there 
are some underlying similarities between these hard X- 
ray emitting symbiotics. Unlike most other symbiotic 
stars, they display a high incidence of optical flickering. 
They also tend to have low optical line strengths, indi- 
cating that they are often only "weakly" symbiotic. Both 
the visibility of the optical flickering (which in most sym- 
bioti cs is overwhelm ed by reprocessed shell-burning emis- 
sion; [Sokoloski 2003) and the weakness or low-ionization- 
state of the optical lines suggest that quasi-steady shell 
burning is not taking place in these objects, either be- 
cause: 1) the WD is more massive, or 2) the accre- 
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tion rate is lower than in other symbiotics. Hard X- 
ray emission might therefore be a proxy for high WD 
mass. In fact, at least one of the hard- X-ray symbiotics, 
T CrB, is a recurrent n ova and contains a high-mass WD 
(|Hachisu fe Katoll2001l ). Finally, jet production appears 
to be more common in flickering symbiotics, and one 
of the 4 hard X-ray symbiotics - CH Cyg - regularly 
produces jets (I Taylor et all 119861 : iKarovska et~aLlll998l : 
ICrocker et afll200U I2002TT Other WDs in hard X-ray 
s ymbiotics might t hus a lso harbor jets. 

ICieslinski et all (|1994D classified RT Cru as a symbi- 
otic star based on its optical spectrum. They noted 
that the lack of strong high-ionization emission lines and 
the very weak forbidden emission lines make the optical 
spectrum similar to that of T CrB. They detected opti- 
cal flickering in the U band with a tim e scale of a few 
tens of minutes. Except for GX 1+4 (jJablonski et al.l 
H997tlChakrabartv fc Roche1H99l . none of the neutron- 
star containing symb i otic s tars produce optical flickering 
fe.g. ISokoloski et all [200lh . which is common in CVs, 
or sho w Balmer or He II emission lines ([Masetti et al.l 
I2007bh . The presence of optical flicker ing and Balmer 
and H e II emission lines from RT Cru (jCieslins ki et alJ 
1994) suggests that it therefore contains an accreting WD 
rather than a neutron star. Reddening estimates from 
optical spectra and infrared magnitudes (coupled with 
the assumptio n that the radius of t he M5 III red giant is 
about 0.5 AU: lvan Belle et al.lll999h suggest that RT Cru 
is between 1.5 and 2 kpc away (J. Mikolajewska, private 
communication) . 

In 2003 and 2004, the IBIS instrument on board IN- 
TEGRAL detected hard X-ray emission extending out 
to ~100 keV from the so urce IGR J12349-6434, which 
IChernvakova et all (l2005h found to hav e a 20- 60 keV flux 
density of ~3 mCrab. iMasetti et alJ (|2005l ) suggested 
an association between IGR J12349-6434 and RT Cru, 
which observations with t he Swift satellite later con- 
firmed ijTueller et al.ll2005t ). The long-term optical light 
curve from the AAVSO indicates that at the time of 
the INTEGRAL observations, RT Cru was in an optical 
bright state; it brightened from 13.5 to 11.5 mag some- 
time between 1998 and 2000. Between 2000 and 2005, 
the optical brightness slowly decreased to approximately 
12.1 mag. The short (4.7 ks) Swift observation of 2005 
August showed that between 2003 and 2005, the hard 
X-ray flux also decreased. 

In this paper, we describe Chandra High Energy 
Transmission Grating (HETG) observations of RT Cru, 
the first member of a new class of hard X-ray emitting 
symbiotic WDs. We detail the observations and data re- 
duction in fj2] and the results from spectral and timing 
analyses in [J3J In 21 we discuss our interpretation of the 
observations, which confirm that the accreting compact 
object in RT Cru is a WD and provide some of the most 
direct evidence to date for an accretion disk around a 
wind-fed WD in a symbiotic system. In this section, we 
also discuss the implications of a system that can accel- 
erate particles to relativistic speeds and produce X-ray 
emission out to greater than 50 keV being powered by 
an accreting WD. We summarize our conclusions in $5] 

2. OBSERVATIONS AND DATA REDUCTION 

On 2005 October 19, the Chandra X-ray Observatory 
performed a 50.1 ks Director's Discretionary Time obser- 



vation of RT Cru using the HETG (|Canizares et al.ll2005h 
and the ACIS-S detector (Obsld 7186, start time 10:21:12 
UT). We requested the DDT observation to attempt to 
catch RT Cru in the optical bright state that appeared 
to be associated with hard X-ray emission. We used the 
HETG instrument because the Swift XRT observation 
of 2005 August hinted at several possible emission-line 
complexes. The data were collected in Timed Exposure 
mode, in which the CCD chips were read out every 2.54 s. 
The data were telemetered back to earth in Faint mode, 
which conveys photon arrival times, event amplitudes, 
and additional information for evaluating the validity 
of each event. We reduced the data according to stan- 
dard procedures using the software package CIAO 3.3 1 . 
We extracted a spectrum from the undispersed light (the 
zeroth-order spot, which fell on the S3, back-illuminated 
chip) using a circular extraction region with a radius of 6" 
centered on the source coordinates: a = 12h 34m 43.74s 
and S = -64° 33'56.0". To obtain the background for the 
zeroth-order light, we extracted photons from a source- 
free circular region on the same CCD. We grouped the 
spectrum, which is shown in Fig. [TJ to have at least 50 
counts per bin. The average zeroth-order source count 
rate was 0.11 c s _1 . 

For the dispersed light from both of the HETG sets of 
gratings - the Medium Energy Grating (MEG) and the 
High Energy Grating (HEG) - we extracted spectra from 
each of the m = ±1, ±2, and ±3 orders individually (us- 
ing the CIAO software tool dmtype2split). To obtain the 
background for the dispersed light, we extracted counts 
from rectangular regions on either side of the spectral 
image. The count rate in the HEG and MEG mil 
orders was 0.042 c s _1 and 0.034 c s _1 , respectively. Al- 
though the dispersed spectral orders contained too few 
counts to produce a high signal-to-noise-ratio spectrum 
of lines spanning the full energy range of the instrument, 
the combined m = ±1 spectrum (grouped at twice the 
full width at half maximum of 0.012A) provided good- 
quality data in the region around the Fe Kq emission-line 
complex. We therefore used the the zeroth-order spec- 
trum for continuum fitting and the dispersed (m = ±1) 
spectra primarily for analysis of the Fe lines. The HEG 
and MEG m = ±2 and ±3 spectra helped confirm the 
Fe line identifications. For spectral fitting of both the 
zeroth and higher order photons, we u sed the standard 
software pack ages Xspec (|Arnaudlll996l ) vl2.3.0 and ISIS 
(Houck 2002). The background contributed less than 1% 
of the total extracted dispersed and undispersed light. 

We generated light curves in the energy bands 0.3- 
4.0 keV and 4.0-8.0 keV by extracting counts (with 
CIAO) from a region containing the zeroth-order spot 
and the m = ±1 dispersed orders of both the HEG and 
MEG. Since we estimated there to be only ~70 back- 
ground counts in this extraction region during the course 
of the observation (compared to more than 9000 source 
counts), we did not background-subtract the light curves. 

At high count rates, two or more photons can arrive 
close enough together in time that they appear to be 
a single event. This "pileup" phenomenon can cause a 
spectrum to become distorted and the count rate to be 
reduced. To confirm that the zeroth order spectrum was 

1 Chandra Interactive Analysis of Observations (CIAO), 
| http : / / exc . harvard . edu/ciao/ 1 
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not affected by pileup, we divided the number of counts 
at the peak of the point spread function of the undis- 
persed light by the number of 2.54-s frames in the obser- 
vation to obtain an upper limit on the number of counts 
per pixel per frame. The resulting 0.08 counts per pixel 
per frame is well below the 1 count per pixel per fram e 
where pileup can become important (|Harris et alJ l2004) . 
Moreover, the average count rate in the zeroth order was 
significantly less than one count per frame time, indicat- 
ing that the light curve was not significantly distorted 
by the saturation that can occur at higher count rates 
(i.e., higher pileup fractions). The pileup fraction in the 
higher-order spectrum was negligible. 

3. ANALYSIS AND RESULTS 
3.1. Spectral Analysis 

To model the X-ray spectrum, we first consider simple, 
single-component continuum models. We fitted these 
models to the binned 0.3 - 8.0 keV zeroth-order spec- 
trum (above 8.0 keV, the noise rises and the quantum 
efficiency drops sharply). Absorbed single-component 
emission models such as a thermal plasma, powerlaw, or 
blackbody (plus Gaussian lines) do not produce accept- 
able fits. Even if we include complex absorption, such as 
an absorber that only partially covers the source, a pow- 
erlaw dist ribution of absorbers (as s een in some magnetic 
CVs; e.g.. lDone fe Magdziarzl[l998h . or a"warm" ionized 
absorber, single-component emission models still do not 
produce acceptable fits. 

Including an additional broad-band emission compo- 
nent improves the fitting results. The spectrum is for- 
mally well fitted with a highly absorbed, optically thin 
thermal plasma (Mekal model in Xspec), plus a moder- 
ately absorbed non-thermal powerlaw component. Since 
there is some degeneracy between the amount of absorp- 
tion and the powerlaw index, we determine the powerlaw 
index by fixing the plasma temperature, T, and fitting 
the spectrum above 4 keV, where absorption is relatively 
unimportant. The resulting photon index is T = 1.05q^, 
where T is given by dF N /dE = K{E/E )- r , dF N /dE 
is the photon flux density, Eq is 1 keV, K is the nor- 
malization constant, and the superscripts and subscripts 
are 90% confidence upper and lower limits, respectively. 
The photon index is not sensitive to the value to which 
we fix T. We determined the remaining model param- 
eters by fixing T to 1.1 and letting the other parame- 
ters vary. The thermal plasma has kT = 8.6l\ 4 keV 
and an absorbing column % = 9.3y 2 2 5 X 10 cm -2 . 
The powerlaw component has an absorbing column of 
Uh = 7.3|'3 x 10 22 cm~ 2 and a normalization K = 
1.2g;g 8 x 10~ 3 photons cm~ 2 s _1 kev" 1 . To obtain an 
acceptable fit to the line emission, the model required 
abundan ces of about 0.3 times so lar (using the abun- 
dances of lAnders fe Grevessdll98"9| ) . Because of the large 
column density absorbing the thermal emission in this 
model, the powerlaw component dominates below ap- 
proximately 4 keV, and the thermal emission dominates 
above 4 keV. We also obtained a formally acceptable 
fit with a highly absorbed powerlaw and a moderately 
absorbed thermal plasma. In that case, the thermal 
plasma dominated the low-energy portion of the spec- 
trum, and the powerlaw dominated the high-energy part 
of the spectrum. 



Finally, we consider the isobaric cooling-flow model 
that has worked well for both magnetic and non-magnetic 
CVs. In this model (Mkcflow; Musho tzkv fe Szvmkowiakl 

1988) , gas is assumed to radiatively cool from a high 
post-shock temperature under conditions of constant 
pressure. Such an isobaric cooling flow has a differen- 
tial emission-measure distribution that is a flat function 
of temperature. The gas is also assumed to be optically 
thin. This cooling flow model provides a good fit to the 
data. We find that the initial post-shock temperature is 
quite high. Although our 0.3-8.0 keV spectrum does not 
allow a high-confidence determination of this parameter, 
the formal 90% confidence lower limit for the initial tem- 
perature is kT max = 55 keV. The minimum cooling-flow 
temperature is consistent with the smallest value allowed 
by the Xspec Mkcflow model (kT m i n = 80.8 eV), indicat- 
ing that the gas does indeed remain optically thin as it 
cools. Allowing the differential emission measure to have 
a powerlaw distribution (with the Cemkl model in Xspec) 
did not improve the fit. As with the two-component con- 
tinuum model, the cooling-flow model requires significant 
absorption. The abundances may al so be sub-solar (0.3 
times solar, using the abundances of lAnders fe Grevessd 

1989) . The best absorption model consists of both a pho- 
toelectric absorber that fully covers the source and an- 
other that only partially covers it. Table [T] lists the best- 
fit parameters for this model, which, as we discuss in fjH 
we believe provides the best description of the Chandra 
spectrum of RT Cru. 

In the first-order spectrum, we detect the iron-line 
complex spanning roughly 6.4 - 7.0 keV. Fig. [5] shows 
the region around the iron-line complex in the combined 
MEG and HEG first-order (to = ±1) spectrum. Because 
of the large absorption and the resulting low count rate at 
low energies, we are not sensitive to lines such as OVIII 
(~ 19 A) and Ne X (~ 12 A) that have been seen in 
HETG observations of some other ac creting WDs (e.g., 
iPandel et all 120051 : iMukai et all 120031) . For the Fe lines, 
we use a simple powerlaw to establish a continuum level 
and three Gaussian profiles to fit the Fe Ka, H-like Fe, 
and He- like Fe lines, respectively. To avoid the possible 
introduction of errors from misalignment of the HEG and 
MEG spectra, we use only the combined HEG first-order 
(to = ±1) spectrum for computation of the equivalent 
widths (EWs). Table [2] lists the line-center energies and 
EWs. Although we do not have enough counts to use the 
recombination, intercombination, and forbidden compo- 
nents of the H- or He-like Fe lines as density diagnostics, 
the observed line strengths and EWs confirm that the 
source is surrounded by a large amount of neutral mate- 
rial and that the abundances might be slightly sub-solar. 
The Fe Ka EW of 108 eV is consistent with that ex- 
pected for a source insid e a cloud of cold material with 
the N ff of - 10 23 cm" 2 (|Inouelll985D . as we found from 
the continuum fitting. 

3.2. Timing Analysis 

Examining time series binned at 508.208 s and 
4065.664 s (i.e., 200 and 1600 times the frame time, re- 
spectively), we detected significant aperiodic, flickering- 
type variations on time scales of minutes to hours in the 
0.3-4.0 keV emission. Figure [3] shows the 508-s binned 
time series (light curves) in the energy ranges 0.3-4.0 keV 
and 4.0-8.0 keV. The fractional amplitude of the stochas- 
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tic variations appears to be largest in the 0.3-4.0 keV 
energy range. In the 508-s binned 0.3-4.0 keV time se- 
ries, the ratio of measured fractional rms variation, s, to 
that expected from Poisson fluctuations alone, s exp , is 
1.96 (s = 36.6% and s exp = 18.7%). We detect the 0.3- 
4.0 keV variability with even greater statistical signifi- 
cance in the 4065-s binned time series; the ratio s/s exp in 
this case is 3.35 (s = 22.1% and s exp = 6.6%). In the 4.0- 
8.0 keV energy range, the 508-s binned time series has 
s/s exp = 1.36 (where s = 16.6% and s exp = 12.2%), and 
the 4065-s binned time series has s/s exp = 1.42 (s = 6.1% 
and s exp = 4.3%). 

We do not detect any periodic flux modulations. We 
are theoretically sensitive to an oscillation with fractional 
amplitude: 

1/2 



expected from accretion, (1/2)(GMM/R) > Lx, the ra- 
dius of the accreting compact object is: 
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where Ctot is the total number of counts in the observa- 
tion (ignoring the small number of background counts, 
which have a negligible effect), 5 and e are small num- 
bers related to the chance that a noise power in the power 
spectrum will exceed the detection threshold (both taken 
to be 0.05), nf r eq is the number of frequencies searched 
(nfreq = 1644), and a has an average value of 0.77 and 
depends upon the locatio n of the signal frequency in 
the fr equency bin (see, e.g.. Ivan der Klisl [i989: Sokoloski 
1999). We are therefore sensitive to oscillations with 
fractional amplitudes of ~ 8% in regions of the power 
spectrum dominated by white noise, which in this case 
consisted of frequencies greater than ~ 1.4 mHz. In this 
analysis, we binned the time series in 15-s bins (i.e., 6 
times the frame time) . We were therefore sensitive to os- 
cillations with periods as short as 30 s, and most sensitive 
to oscillations with periods between 30 s and 12 m. 

Consistent with the presence of flickering in the light 
curves, the power spectrum rises at frequencies below 1.4 
mHz. At these low frequencies, the power spectrum has 
a powerlaw index (i.e., slope on a log- log plot) of about 
— 1. This "l//-noise" at low frequencies reduces our sen- 
sitivity to oscillatory signals with periods greater than 
approximately 12 m. Since we expect the minimum os- 
cillation amplitude to which we arc sensitive to increase 
roughly as the square root of the rising average broad- 
band power as we go to lower frequencies (i.e., longer 
periods), the oscillation amplitude required for detection 
increases gradually from >8% to >15% as we move from 
periods of 12 m to 1 hr. Taking into account the underly- 
ing broadband power, as well as the number of frequency 
bins searched, we did not detect any statistically signifi- 
cant oscillations in any portion of the power spectrum. 

4. DISCUSSION 

4.1. Interpretation of the Chandra Observations 

To estimate the radius of the accreting compact object, 
we take the unabsorbed 0.3-8.0 keV X-ray luminosity, 
Lx, to be either approximately equal to, or a rough lower 
limit to, the emission from an accretion-disk boundary 
layer (we justify the assumption that the Chandra X- 
ray emission emanates from a boundary layer in the para- 
graphs that follow). Comparing this luminosity with that 



R < 3.2 x 10 8 cm 




M 



1.8 x 10- 9 M Q yr- 



( 3 ) 

where R and M are the radius and mass of the accre- 
tor, respectively, and M is the rate of accretion through 
the boundary layer. The radius is that of a WD. The 
Chandra X-ray spectrum therefore confirms that the 
compact object is a WD. 

To determine whether the Chandra-b&iid X-ray emis- 
sion is indeed from an accretion-disk boundary layer, we 
consider the rapid variability. Rapid flickering typically 
emanates from an accretion region close to a compact ob- 
ject. Our detection of nickering therefore suggests that 
the X-ray emission detected from RT Cru by Chandra is 
powered by accretion. This accretion could proceed via a 
wind-fed accretion disk, magnetic accretion columns, or 
Bondi-Hoyle type direct-impact of the accreting material 
onto the WD. While the two-component (thermal plasma 
plus powerlaw) model provides a formally acceptable fit 
to the data, it is difficult to construct an interpretation 
of this model that is consistent with the rapid flicker- 
ing from accretion onto a WD. The isobaric cooling-flow 
spectral model, on the other hand, 1) provides a good fit 
to the data, 2) has been successfully applied to both the 
boundary layer emission from non-m agnetic CVs (e.g., 
iMukai et al]|2003t iPandel et al"1l2005h and th e accretion 
columns of magnetic CVs fe.g. JCropper et al.lll998l ) , and 
3) provides a natural context for the flickering from ac- 
cretion. 

Most CVs with X-ray emission as hard as that 
which INTEGRAL and Swift /BAT have detected from 
RT Cru have magnetic fields strong enough to channel 
the accretion flow into accretion columns (B ~ 10 5 - 
10 6 G, where B is the magnetic field strength at the 
WD surface). Of the 8 CVs detected at energies greater 
than - 50 keV with Swi ft /BAT, all but S S Cyg are 
likely magnetic accretors (Barlow et al.l 12006). In these 
systems, the hard X-rays come from hot gas behind the 
stand-off shock in the accretion column. Magnetic CVs 
typically have X-ray oscillations with pulsation ampli- 
tudes of tens of percent at t he WD spin p eriod, which 
is usually less than an hour (| Warned 1 1995h . Since typi- 
cal symbiotic-star accretion rates are higher than typical 
CV accretion rates, if the WD in RT Cru was strongly 
magnetic and in spin equilibrium, the spin period would 
probably be either comparable to or faster than those 
in CVs. Given our sensitivity to oscillations with peri- 
ods less than an hour, we therefore should have detected 
a spin modulation if RT Cru was magnetic. In fact, the 
power spectrum has no statistically significant peaks. We 
conclude that RT Cru is probably not a magnetic accre- 
tor. We thus favor the picture in which the X-ray emis- 
sion from RT Cru detected by Chandra is from a cooling 
flow in an accretion-disk boundary layer. 

4.2. Implications 

The parameters of the cooling-flow fit to the boundary- 
layer emission provide the accretion rate as well as infor- 
mation about the WD. From the normalization param- 
eter of the cooling-flow model, the accretion rate onto 
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the WD is M = 1.8 x 1(T 9 M yr" 1 (d/2kpc) 2 (see Ta- 
ble Q}. Since quasi-steady nuclear burning is frequently 
present on the WD surface, symbiotic stars should have 
accretion rates that are on average higher than those in 
CVs. The accretion rate we have found for RT Cru is 
consistent with this picture. It is also, however, just low 
enough that we expect the bound a ry lay er to remain op- 
tically thin; iNaravan fe Pophaml (|1993f ) find that for a 
1 Mq WD, the boundary layer remains optically thin for 
accretion rates below 3 x 10~ 9 M© yr _1 . 

The parameters of the cooling-flow fit also indicate 
that the WD radius is small, suggesting that the WD 
could be quite massive. Taken at face value, the ra- 
dius constraint would imply that the WD mass is at 
least 1.3 Mq. The high upper cooling-flow tempera- 
ture, kT max > 55 keV (see TableQ]), supports the conclu- 
sion that the accretor is a massive WD. The relationship 
between kT max and WD mass is due to the fact that 
that the Kepler velocity (v^ = GM/R) is greater in the 
deep potential well of a more massive WD. Since the 
boundary-layer material is shock heated, and the initial 
post-shock temperature (T max ) is proportional to veloc- 
ity squared, T max increases with WD mass. Alterna- 
tively, if we equate the amount of energy available per 
particle, (l/2)/im p w^, where is the mean molecular 
weight and m p is the mass of a proton, with the energy re- 
leased per particle in an isobaric cooling flow, (5/2)kT max 
()Pandel et al.ll2005[ ). we see that kT max oc u| oc GM/R. 
Although the determination of kT rnax from X-ray emis- 
sion below 8 keV is highly uncertain, we can still ask 
what such a high kT max would imply if it is confirmed 
by an instrument with greater high-energy sensitivity 
(such as SUZAKU) . For t heir sample of 9 non-magnetic 
CVs, iPandel et all ((III) found the upper cooling-flow 
temperature was roughly consistent with the expected 
kT max = (3/5) kT V i r , where T V i r is the virial temperature 
(defined by (3/2)kT vir = ( l/2)um p vl). Using the W D 
mass-radius relationship oLHa nsen fc Kawalerl (|1994f ). a 
maximum cooling-flow temperature of kT max > 55 keV 
implies M > 1.3 M . 

The high absorbing columns for both the partial- 
covering and fully covering absorber, as well as the cover- 
ing factor of >0.7 for the partial-covering absorber (see 
Table Q}, indicate that the X-ray source is highly ob- 
scured at this epoch. Since the X-ray emission region 
is small, the absorber that only partially covers the X- 
ray source must also be small. A possible source of this 
partially covering absorber is an accretion structure such 
as an accretion dis k seen almost edge-on (as in OY Car; 
IPandel et all [20051. We assume that the fully covering 
component of the absorption comprises both interstel- 
lar absorption (1.1 x 10 22 cnr 2 from NASA/IPAC IRSA) 
and intrinsic absorption. The column density of this ab- 
sorber is probably high because the WD orbits within the 
strong, dense stellar wind from the red giant. RT Cru 
has an orbital period of ~ 450 d (J. Mikolajewska, pri- 
vate communication) , and therefore a binary separation 
on the order of an AU. This separation puts the WD well 
within the dense region of the red-giant wind. Month- 
time-scale variations in the absorption (not correlated 
with the orbital pe riod) have been detected by Swift 
(jKennea et al.ll2007l ). suggesting that either the red-giant 
wind is clumpy, that the mass loss rate in the wind of the 
red giant is variable, or the accretion structure or struc- 



tures partially blocking the WD boundary layer must be 
unstable. 

4.3. The Nature of RT Cru 

As shown by the INTEGRAL detection in 2003-2004 
(jChernvakova et al.ll2005l ). RT Cru can at times produce 
X-ray emission out to greater than 60 keV. From the 
broad -band flux densities reported bv lChernvakova et al.l 
(2005), the hard X-ray spectrum in 2003 - 2004 appears 
consistent with a power law with photon index T = 2.7 
(energy index a — 1.7), and the 16-100 keV luminos- 
ity at that time was approximately 10 Lq (d/2kpc) 2 . By 
2005, however, the 16-100 keV luminosity had dropped, 
and the hard X-ra y spectrum was closer to thermal 
(|Kennea et aLll2007t ). To better appreciate the proper- 
ties of the unusual accreting WD in RT Cru, we briefly 
explore the possible sources of the powerlaw hard X-ray 
emission observed by INTEGRAL in 2003-2004. In par- 
ticular, we consider direct synchrotron emission, inverse- 
Compton (IC) scattering from a thermal distribution of 
electrons, and IC scattering from a non-thermal distri- 
bution of electrons. IC scattering from a non-thermal 
distribution of electrons turns out to be the most likely 
option. 

One way to generate a powerlaw energy spectrum 
is with synchrotron emission from a powerlaw distri- 
bution of relativistic electrons moving in a magnetic 
field. Most of the synchrotron emission from an elec- 
tron with Lorentz factor 7 is emitted at a frequency 
v S ync = (0.3/27r)(3sina/2)7 2 ync (gi3/m e c), where a is 
the pitch angle, q is the electron charge, B is the mag- 
netic field stre ngth, m e is the electron ma ss, and c is the 
speed of light (|Rvbicki fc Lightmanl lT979). Since we did 
not detect X-ray pulsations from RT Cru, the magnetic 
field at the surface of the WD is probably not strong 
enough to disrupt the accretion flow, and therefore less 
than ~ 10 4 G. Thus, even near the surface of the WD, 
Lorentz factors r y syn c of a few times 10 4 to 10 5 would be 
needed to produce significant direct synchrotron emission 
at 60 keV. We can estimate the maximum Lorentz fac- 
tor to which electrons will be accelerated by setting the 
diffusive shock acceleration time scale equal to the syn- 
chrotron cooling tim e scale. Following the approach of 
iMarkoff et al.1 (|2001l ). we equate the syn chrotron loss r ate 
to a conservative acceleration rate (see lJokipiilll987t ) to 
find a maximum Lorentz factor to which electrons can be 
accelerated of lmax w 2700(B/10 4 G)- 1 / 2 (C/100)- 1 / 2 , 
where £ is the ratio of the diffusive scattering mean free 
path to the gyroradius. Thus jmax is one to two orders 
of magnitude below ^sync- 

Moreover, whereas significant 
power from a distribution of synchrotron-emitting rela- 
tivistic electrons would also be expected at energies be- 
low 15 keV, the 16-100 keV X-ray luminosity of ~ 10 Lq 
during 2003 and 2004 is already close to the total en- 
ergy budget available from accretion. Finally, to gen- 
erate the X-ray spectral index observed by INTEGRAL 
from direct synchrotron emission, one would need a dis- 
tribution of electrons that is much steeper than expected 
from standard theor ies of particle acceleration in shocks 
(|Ellison et al.l [19901 ). It is therefore unlikely that the 
hard X-ray emission from RT Cru was due to direct syn- 
chrotron emission. 

As an aside, we note that the production of significant 
synchrotron emission at radio wavelengths would not re- 
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quire such high Lorentz factors. The ratio of synchrotron 
power to IC power from an electron is equal to the ratio of 
the energy density in the magnetic field to energy density 
in photons (|Rvbicki fc Lightmanlll979f ). Examining this 
ratio as a function of distance from the WD, we expect a 
similar amount of power from direct synchrotron and IC 
scattering near the surface of the WD, where the B field 
could be strong (B ~ 10 4 G). We therefore suggest that 
the next time RT Cru is in a powerlaw hard X-ray state 
like the one detected by INTEGRAL in 2003-2004, radio 
observations be performed to look for radio synchrotron 
emission. 

Hard X-rays can be produced with much lower Lorentz 
factors through IC scattering. For p hotons scattering off 
of a th ermal distribution of electrons. iRevnolds fc Nowald 
(2003) give a relation between the observed photon in- 
dex r and the Compton y parameter, which is related 
to the factor by which the average photon energy in- 
creases. A steep photon index of T > 2 indicates a y 
parameter less than 1, in which case there is no sig- 
nificant up-scattering. Therefore, since RT Cru had a 
steep photon index of T w 2.7, scattering off of a ther- 
mal distribution of electrons could not have been re- 
sponsible for the powerlaw hard X-ray emission. Non- 
thermal, relativistic electrons must have been involved 
in producing the powerlaw spectrum observed by INTE- 
GRAL. RT Cru thus contains a white dwarf that can 
accelerate electrons to relativistic speeds. Three other 
systems that contain WDs that also generate relativis- 
tic electrons are CH C yg, RS Oph, and possibly R Aqr, 
all of which have jets (iRupen et a,l.ll2007t lO'Brien et al.l 
[20061 iCrocker et al.ll2001i iNichols et al.ll2007h . 

Assuming now that the powerlaw hard X-ray spectrum 
was due to IC scattering from a powerlaw distribution of 
relativistic electrons, we can estimate the location of the 
scattering electrons. Since the strength of a dipole B field 
falls like 1/r 3 , whereas the energy density in the photon 
field only falls like 1/r 2 , IC scattering will dominate as 
one moves farther from the WD. If we ask how far away 
from the radiation source the IC region should be to give 
us an IC cooling time on the order of a year (the approx- 
imate duration of the powerlaw hard X-ray state), we 
find that it should be a few tenths of an AU away. Given 
the orbital period for this systems, that could put the 
IC emission region either between the WD and the red 
giant, as in a colliding-winds region, above the WD disk 
in a corona, or at the base of a jet. A model consisting 
of IC scattering off of relativistic electrons in a corona 
or at the base of a jet reproduces the observed hard X - 
ray emission well in X-ray binaries (|Markoff et al.ll2005h . 
The steep powerlaw index, which is similar to that seen in 
the steep powerlaw state in microq uasars, could be due 
to a low scattering optical depth (jRvbicki fc Lightmanl 
[19791) . 



5. CONCLUSIONS 

We have observed the first symbiotic star with X-ray 
emission out to greater than 60 keV with the HETG 
on Chandra. The stochastic variability and cooling- flow 
type spectrum suggest that RT Cru is powered by ac- 
cretion onto a WD through a disk with an optically thin 
boundary layer. The accretion rate is near the top of the 
range for which the boundary layer can remain optically 
thin. The high initial temperature of the cooling flow, 
and the high luminosity given the accretion rate from the 
spectral fit suggest that the WD in RT Cru could be quite 
massive. More generally speaking, it would be difficult to 
get such hard X-ray emission from accretion onto a low- 
mass WD with a shallow potential well and low Kepler 
velocity. Given the nature of the powerlaw hard X-ray 
emission previously observed by INTEGRAL, it therefore 
appears that the accreting, non-magnetic WD in RT Cru 
is able to generate a non-thermal, powerlaw distribution 
of electrons and very hard X-ray emission through IC 
scattering. Three other systems in which WDs can accel- 
erate electrons to rela tivistic speeds (RS Oph, CH Cyg, 
and possibly R Aqr: IRupen et al.l 120071 : lO'Brien et alJ 
[20061 ICrocker et all 120011: INichols et al.l 120071) all have 
jets. Radio observations of RT Cru during the next pow- 
erlaw hard X-ray state like that observed by INTEGRAL 
in 2003 and 2004 could play an important role in diag- 
nosing this system, and determining the extent to which 
some symbiotic stars m ight constitute the nanoqu asar 
analog to microquasars (jZamanov fc Marzianill2002n . 
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TABLE 1 
Cooling-flow model parameters. 



Parameter Value (Min, Max) 1 

M b (It)" 9 M yr" 1 ) 1.8 (1.6, 2.0) 

kT max (keV) 80 (56, ■ ■ ■ ) 

n H : Full covering (10 22 cm" 2 ) 8.2 (7.7, 8.8) 

n H : Partial covering (10 22 cm' 2 ) 65 (52, 78) 

Covering Fraction 0.74 (0.69, 0.79) 

Abundance (w.r.t. Solar c ) 0.30 (0.02, 0.45) 

F x d (10~ 12 erg cm" 2 s" 1 ) 9.1 (7.5, 10.2) 

L x d (10 34 er S s_1 ) 3.1 (2.4, 3.9) 

Note. — The model consists of optically thin thermal emission from an isobaric cooling flow with absorbers that both fully cover and 
partially cover the source, plus a Gaussian line. 

a 90% confidence upper and lower limits. b Accretion rate onto the compact object. Using solar abundance of lA"ndcrs & G revessel 
Ill989l ). d F x is the absorbed 0.3-8.0 keV flux, and L x is the unabsorbed 0.3 - 8.0 keV luminosity (d = 2 kpc). 
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TABLE 2 
Iron lines. 





Fe XXV 


Fc XXVI 


Fe Ka 


Line center (keV) a : 
EW b (eV): 


6.946| : !i 

72 


6.693|;™6 
60 


6.379| : i! 
108 



a Super- and subscripts represent 90% confidence upper and lower limits, respectively. 15 Gaussian-fit equivalent widths. EW uncertainties 
arc on the order of 10-15%. 
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Energy (keV) 

Fig. 1. — Undispersed (zeroth order) spectrum. The top panel shows the spectrum with the absorbed, isobaric cooling-flow model 
over-plotted. The bottom panel shows residuals with respect to this model (in units of x 2 , where x 2 is shorthand for the difference between 
the data and the model, squared, divided by the variance, with the sign of the difference between the data and the model). 
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Fig. 2. — The iron-line complex from the combined HEG and MEG first-order (m = ±1) spectrum. The best fit model of a powerlaw 
plus three Gaussian emission lines is over-plotted. The bottom panel shows the residuals, in the same units as Fig. [T] 
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Fig. 3. — Chandra light curves for RT Cru, with a bin size of 508.28 s. The light curves include the undispersed light as well as the 
counts from the HEG and MEG m = ±1 orders. The top and bottom panels show the flux as a function of time in the energy ranges 
0.3—4.0 keV and 4.0-8.0 keV, respectively. The 0.3—4.0 keV emission is clearly variable on time scales of minutes to hours. 



